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ABSTRACT 

Galaxy clustering data can be used to measure the cosmic expansion history H(z), the 
angular-diameter distance Da(z), and the linear redshift-space distortion parameter /?(z). 
Here we present a method for using effective multipoles of the galaxy two-point correlation 
function (£o(s)> &( s )> £,4(s), and £e(s)> with s denoting the comoving separation) to measure 
H(z), Da(z), and and validate it using LasDamas mock galaxy catalogs. Our definition 
of effective multipoles explicitly incorporates the discreteness of measurements, and treats the 
measured correlation function and its theoretical model on the same footing. We find that for 
the mock data, £o + £2 + £4 captures nearly all the information, and gives significantly stronger 
constraints on H(z), Da(z), and f3(z), compared to using only £0 + £2- 

We apply our method to the sample of luminous red galaxies (LRGs) from the Sloan Dig- 
ital Sky Survey (SDSS) Data Release 7 (DR7) without assuming a dark energy model or a flat 
Universe. We find that £i(s) deviates on scales of s < 60 Mpc/h from the measurement from 
mock data (in contrast to £o(s)> £2(3), and £e(s))> leading to a significant difference in the 
measured mean values of H(z), Da(z), and f3(z) from £0 +£2 and £0 + £2 +£4, thus it should 
not be used in deriving parameter constraints. We obtain {H(0. 35), -0^(0.35), n, m h 2 , f3(z)} 
= {Tg.et^kms^Mpc- 1 , 1057±^Mpc, 0.103 ± 0.015, 0.44 ± 0.15} using | + | 2 - 
We find that H(0. 35) r s (zd)/c and DA(0.35)/r s (zd) (where r s (zd) is the sound horizon 
at the drag epoch) are more tightly constrained: |iJ(0.35) r s (zd)/c, D y i(0.35)/r s (z ( ;)} = 

{0.0437;°;oo43,6.48;°;43} using £ + 6- We conclude that the multipole method can be used 
to isolate systematic uncertainties in the data, and provide a useful cross-check of parameter 
measurements from the full correlation function. 
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1 INTRODUCTION 

The cosmic large-scale structure from galaxy redshift surveys pro- 
vides a powerful probe of dark energy and the cosmological model 
that is h ighly complementary to the cosmic micr owave backg round 
(CMB) ^Bennett et al|2003h . supernovae (S Ne) teiess et alj | l998j : 
Perlmutter et al.l 1 19991) . and weak lensing ( Wittman et al.i |200I 



Bacon. Refregier. & EUisl|200Ct Ik aiser. Wilson. & Luppinol 



2000; 



Van Waerbeke et al.ll2000r) . The scope of galaxy redshift surveys 



has dramatically increased in the last decade. The PSCz surveyed ~ 
15, 000 galaxies usin g the Infrared Astronomical Satellite (IRAS) 
dSaunders et atfeoOOh . the 2dF Gal axy Redshift Survey (2dFGRS) 
obtained 221,414 galaxy redshifts dColless et al1l200Ul2003l) , and 



the Sloan Digital Sky Survey (SDSS) has c ollected 930,000 galax y 
spectra in the Seventh Data Release (DR7) lAbazaiian et alj2009 ) . 
The ongoing galaxy surveys will probe the Universe at higher 
redshifts; WiggleZ is surveying 240,000 e mission-line galax ies at 
0.5 < z < 1 over 1000 square degrees dBlake et al.ll2009t) . and 
BOSS is surveying 1.5 million luminous red galaxies (LRGs) a t 
0.1 < z < 0.7 over 10,000 square degrees (Eise nstein et al.l201ll) . 
The planned space mission Euclid^ will survey over 60 million 
emission-line galax i es at 0.7 < z < 2 over 15,000 square degrees 
dCimatti et al.ll2009l ; lwa"ng et al J20 id : iLaureii s et al.ll201 lh . 

Large-scale structure data from galaxy surveys can be an- 
alyzed using either the power spectrum or the correlation func- 
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tion. Although these two methods are simple Fourier transforms 
of one another, the analysis processes are quite different and 
the results cannot be converted using Fourier transform directly 
because of the finite size of the survey volume. The SDSS 
data have been analyzed using both the power spectrum method 
(see, e.g.. rTegmark et alj |2004 iHutsil 120051; Padmanabhan et al 



selection criteria were used in previous paper s ( Zehavi et alJI 20051 : 
lEisenstein eTail 120051; lOkumura et all 120081: iKazin et alj l2010ah. 



20071 : iBlake et alj 120071 : IPercival et all 120071 |2010| ; iReid et al 
2010; Montesa noetaiT 1201 ll ). " and the correlation function 



method (see, e.g., lEisenstein et alj 120051; lOkumura et alj 120081; 



Cabre & Gaztanaga 2009; Martinez et al. 2009; Sanche z et all 
20091: IKazin et alj |20TM; IChuang, Wang. & Hemanthj 1 201% 



Samushia etaTlb oi 1 ; Padmanabhan et al. 20121) . 



The power of galaxy clustering as a dark energy probe lies 
in the fact that the Hubble parameter, H(z), and the angular di- 
ameter distance, Da(z), can in principle be extracted simultane- 
ously from data through the measurement of the baryon acous- 
tic oscillation (BAO) scale in the radial and transverse direction s 
dBlake & Glazebrook|l200l ISeo & Eisens"teinll2003l ; IWang||2006t) . 
lOkumura et al.l d2008l) concluded that SDSS DR3 LRG data were 
not sufficient for measuring H(z) and Da{z); they derived con- 
straints on cosmological parameters assuming that dark energy is 
a cosmological constant. ICabre & Gaztanagal ([2009 ) measured the 
linear redshift space di stortion parameter P, galaxy bia s, and erg 
from SDSS DR6 LRGs. lGaztanaga. Cabre. & Huild2009}) obtained 
a measurement of H(z) b y measuring the pea k of the 2PCF along 
the line of sight. However. lKazin et alj (2010b) showed that the am- 
plitude of the line-of-sight p eak is consistent with sa mple variance. 

In our previous paper dChuang & WandlioT l|), we presented 
the method to measure H(z) and Da(z) from t he full 2D corre- 
lation function of a sample of SDSS DR7 LRGs l lEisenstein et alj 
120011) without assuming a dark energy model or a flat Universe. It 
is also the first application which inclu des the geometric distortio n 
(also known as Alcock-Paczynski test dAlcock & Paczvnskll979h ) 
on the galaxy clustering data at large scales. We demonstrated the 
feasibility of extracting H(z) and Da(z) by applying the method 
on the individual LasDamas mock catalogs which mimic the galaxy 
sample and survey geometry of the observational data we used. In 
this paper, we extend our method by exploring the use of the mul- 
tipoles of the correlation function to measure H(z), Da(z), and 
/3(z). The obvious advantage of using multipoles of the correlation 
function instead of the full 2D correlation function is the reduced 
number of data points used to obtain similar amount of information. 
In Section[2] we introduce the galaxy sample used in our study. In 
Section [3] we describe the details of our method. In Section [4] we 
present our results. In Sectioin[5] we apply some systematic tests to 
our measurements. We summarize and conclude in Sec. [6] 



2 DATA 

The SDSS has observed one-quarter of the entire sky and per- 
formed a redshift survey of galaxies, quasars an d stars in five 
passbands u, q, r. i, and z w ith a 2.5m telescope dFukugita et alj 
ll996l : lGunn et al.ll998tl2006h . We use th e public catalog, the NYU 
Value-Added Galaxy Catalog (VAGC) dBlanton et afll2005ri . de- 
rived fr om the SDSS II final public data release, Data Release 7 
(DR7) dAbazaiian et We select our LRG sample from 

the NYU VAGC with the flag primTarget bit mask set to 32. 
K-corrections have been applied to the galaxies with a fiducial 
model (ACDM with f2 m = 0.3 and h = 1), and the selected 
galaxies are required to have rest-frame q-band absol ute magni- 
tudes -23.2 < M g < -21.2 dBlanton & Roweij2007h . The same 



The sa mple we use is referred to as "DR7full" in Kaz in et all 
d2010id) . Our sample includes 87000 LRGs in the redshift range 
0.16-0.44. 

Spectra cannot be obtained for objects closer than 55 arcsec 
within a single spectroscopic tile due to the finite size of the fibers. 
To correct for these "collisions", the redshift of an object that failed 
to be measured would be assigned to be the same as the nearest 
successfully observed one. Both fiber collisi on corrections and K - 
corrections have been made in NYU- VAGC (Blanto n et al. 2005). 
The collision corr ections applied here are different from what has 
been suggested in IZehavi et alj d2005l) . However, the effect should 
be small since we are using relatively large scale which are less 
affected by the collision corrections. 

We construct the radial selection function as a cubic spline fit 
to the observed number density histogram with the width Az = 
0.01. The NYU- VAGC provides the description of the geometry 
and completeness of the survey in terms of spherical polygons. We 
adopt it as the angular selection function of our sample. We drop 
the reg ions with complet eness below 60% to avoid unobserved 
plates dZehavietai]|2005h . The Southern Galactic Cap region is 
also dropped. 



3 METHODOLOGY 

In this section, we describe the measurement of the multipoles of 
the correlation function from the observational data, construction 
of the theoretical prediction, and the likelihood analysis that leads 
to constraints on dark energy and cosmological parameters. 



3.1 Measuring the Two-Dimensional Two-Point Correlation 
Function 

We convert the measured redshifts of galaxies to comoving dis- 
tances by assuming a fiducial model, ACDM with Q m — 0.25. 
We use the two-point correlation function estimator given by 
lLandv&Szalavldl993l) : 



DD(a, tt) - 2DR(a, ir) + RR(a, n) 
RR(a, tt ) 



(1) 



where n is the separation along the line of sight (LOS), a is the sep- 
aration in the plane of the sky, DD, DR, and RR represent the nor- 
malized data-data, data-random, and random-random pair counts 
respectively in a distance range. The LOS is defined as the direc- 
tion from the observer to the center of a pair. The bin size we use 
here is 1 /i _1 Mpcx 1 /i _1 Mpc. The Landy and Szalay estimator has 
minimal variance for a Poisson process. Random data are generated 
with the same radial and angular selection functions as the real data. 
One can reduce the shot noise due to random data by increasing the 
number of random data. The number of random data we use is 10 
times that of the real data. While calculating the pair counts, we as- 
sign to each data point a radial weight of 1/[1 + n[z) ■ P w ], where 
n(z) is the radial selec tion function and P w — 4 ■ 10 4 /i~ 3 Mpc 3 
dEisenstein et alj2005l) . 

3.2 Theoretical Two-Dimensional Two-Point Correlation 
Function 

We compute the linear power spectra by using CAMB 
dLewis . Challino r. & Lasenb v hoooh . To include the effect of non- 
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linear structure formation on the BAOs, we first calculate the 
dewiggled power spectrum 



Pdw(k) = P iin (fc)exp ( -7^2 )+Pnm(k) 



1 



exp 



~2k 2 



where Pu„ (k) is the linear matter power spectrum, P nw (k) is the 
no-wi ggle or pure CDM pow er spectrum calculated using Eq.(29) 
from lEisenstein & Hul dl998l) . and fc* is marginalized over with a 

flat prior over the range of 0.09 to 0. 13. 

We then use the software package halofit dSmith et alJliooH) 
to compute the non-linear matter power spectrum: 



spectra. On the scales we use for comparison with data, the theoreti- 
cal correlation function only depends on cosmic curvature and dark 
energy through parameters H(z), Da(z), and P(z), assuming that 
, (2) dark energy perturbations are unimportant (valid in simplest dark 
energy models). Thus we are able to extract constraints from data 
that are independent of a dark energy model and cosmic curvature. 



rhaiofit{k) 
Pm{k) 



* halo fit, riw (A;) 

(k) 

Pdw(k)r ha iofit{k), 



(3) 
(4) 



where Phalofu,nw{k) is the power spectrum obtained by applying 
halofit to the no-wiggle power spectrum, and P n i(k) is the non- 
linear power spectrum. We compute the theoretical real space two- 
point correlation function, £(»"), by Fourier transforming the non- 
linear power spectrum P n i(k). 

In the linear regime (i.e., large scales) and adopting the small- 
angle approximation (which is valid on scales of interest), the 2D 
corre l ation function in the redshift space can be written as dKaisej 
ll987l:lHamiltor]|l992h 



(5) 



where s — \/a 2 + tv 2 , fi is the cosine of the angle between 
s = (a, tt) and the LOS, and Pi are Legendre polynomials. The 
multipoles of £ are defined as 



Co(r) = 
6(r) = 



, 2/3 p 2 
1 + — + — 

3 5 



f + ^W)-f(0], 



Hr) + ~£» - ~f (r) 



Ur) = -35 
where f3 is the redshift space distortion parameter and 



(6) 
(7) 
(8) 

(9) 
(10) 



Next, we convolve the 2D correlation function with the distribution 
function of ran dom pairwis e velocities, f(v), to obtain the final 
model £0, tt) dPeeblesll 1980h 



£(r) = ^ i{r')r' 2 dr\ 



£(r) = — £,{r')r' 4 dr' 



£(<7,7r) = / £ cx,7r 



H(z)a(z) 



f(v)dv, 



(11) 



where the random motions are rep resented by an exponential form 
dRatcliffe et al.ll 19981 : lLandvll2002h 



/(«) 



a v \/2 



exp 



V2\v\ 



(12) 



where a v is the pairwise peculiar velocity dispersion. 

The parameter set we use to compute the theoretical cor- 
relation function is {H(z), Da(z), ft, Q m h 2 , ftbh 2 , n s , a v , fc*}, 
where Q. m and il;, are the density fractions of matter and baryons, 
n s is the powerlaw index of the primordial matter power spec- 
trum, and h is the dimensionless Hubble constant (Ho ~ lOO/i km 
s _1 Mpc _1 ). We set h = 0.7 while calculating the non-linear power 



3.3 Effective Multipoles of the Correlation Function 

From Eqs.(f5} and dl It . we find that 



dvf(v)i 




H(z)a(z) 



21 + 1 
2 

21 + 1 



dfJ.£(cr y ir)Pi(fi) 



dBy/l-lJ?Z(*,ir)Pi(n), 



(13) 



where fi = cos 9, and Pi (fj,) is the Legendre Polynomial (I =0, 2, 
4, and 6 here). Note that we are integrating over a spherical shell 
with radius s, while actual measurements of £(cr, 71") are done in 
discrete bins. To compare the measured £(er, n) and its theoretical 
model on the same footing, we convert the last integral in Eo.lll3t 
into a sum. This leads to our definition for the effective multipoles 
of the correlation function: 



E 



(2/ + l)C(a,7r)P i (M)v /l -M 2 



^ ^ ^ Number of bins used in the numerator ' ^ 

where As = 5 /i _1 Mpc in this work, and 



a = ( n + 2^ h ~ lM P c ' "- = 0,1, 2, ... 
7r = (m + i)/i _1 Mpc, m = 0, 1, 2, . 



/'■ : - 



V<7 2 + TT 2 



(15) 
(16) 
(17) 



Note that both the measurements and the theoretical predic- 
tions for the effective multipoles are computed using Eq. (114l >. with 
f (a, 71") given by the measured correlation function (see Eq.([T}) for 
the measured effective multipoles, and Eqs.(f5}-([TTJ for their the- 
oretical predictions. We do not use the conventional definitions of 
multipoles to extract parameter constraints as they use continous 
integrals. Bias could be introduced if the definitions of multipoles 
are diffirent between measurements from data and the theoretical 
model. 



3.4 Covariance Matrix 

We use the 160 mock catalogs from the LasDamas simulations^ 
(McBride et al., in preparation) to estimate the covariance matrix 
of the observed correlation function. LasDamas provides mock cat- 
alogs matching SDSS main galaxy and LRG samples. We use the 
LRG mock catalogs from the LasDamas gamma release with the 
same cuts as the SDSS LRG DR7full sample, -23.2 < M g < 
—21.2 and 0.16 < z < 0.44. We have diluted the mock catalogs to 
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match the radial selection function of the observational data by ran- 
domly selecting the mock galaxies according to the number density 
of the data sample. We calculate the multipoles of the correlation 
functions of the mock catalogs and construct the covariance matrix 
as 



N 



IS 



(18) 



where N is the number of the mock catalogs, X m is the mean of 
the m th element of the vector from the mock catalog multipoles, 
and X^ is the value in the m th elements of the vector from the k th 
mock catalog multipoles. The data vector X is defined by 



l SO ! S2 > S2 i 



i S2 ' 



(19) 



where N is the number of data points in each measured multipole; 
N — 16 in this work. The length of the data vector X depends on 
how many multipoles are used. 



3.5 Likelihood 

The likelihood is taken to be proportional to exp(— X 2 /2) 
jPress et aill992l) , with \ 2 given by 




s (Mpc/h) 



Figure 1. Measurement of monopole of the correlation function of the 
SDSS DR7 LRG (diamond data points) comparing to the averaged 
monopole of the correlation functions of the mock catalogs (solid line). 
The error bars are taken as the square roots of the diagnal elements of the 
covariance matrix. 



^ ] [Xth,i ~ Xobs.i] Cij [Xth.j — X i 



bs,j\ 



(20) 



* j =1 



where Nx is the length of the vector used, Xth is the vector from 
the theoretical model, and X a b s is the vector from the observational 

data. 

As explained in IChuang&WanglfcOllI) . instead of recalculat- 
ing the observed correlation function for different theoretical mod- 
els, we rescale the theoretical correlation function to avoid render- 
ing x 2 values arbitrary. The rescaled theoretical correlation func- 
tion is computed by 

where fth is given by eq. i I 1 [ i. Hence \ 2 can be rewrited as 

x = { T 1 Xth,* - xli*^ } c ft l d ij ■ 



{T^X, 



th.J 



X 



fid 

obs.j 



} 



(22) 



where T~ 1 Xth is a vector given by eq. i !2 1 1 with £th replaced 
by its effective multipoles (defined by eq. d!4t). and X^ is the 
corresponding vector from observational data measured assum- 
ing the fiducial model in converting redshifts to distances. See 
IChuang & Wansl feoilh for a more detailed description of our 
rescaling method. 



3.6 Markov Chain Monte-Carlo Likelihood Analysis 

We use CosmoMC in a Mark ov Chain Monte-Carlo likelihood 
analysis jLewis & Bridiel [2002). The parameter space that we ex- 
plore spans the parameter set of {H(0.35), Da (0.35), 0, m h 2 , /3, 
n b h 2 , n„, er v , fc*}. Only {#(0.35), £>a(0.35), n m h 2 , /3} are 
well constrained using SDSS LRGs alone in the scale range in- 
terested. We marginalize over the other parameters, {Qth 2 , n s , 
ct v , k*}, with the flat priors, {(0.01859,0.02657), (0.865,1.059), 
(0,500)s _1 km, (0.09,0.13)/iMpc" 1 }, where the flat priors of 



Qbh and n s are centered on the measurements from WMAP7 
and ha s width of ±7ctwmap (with owmap from lKomatsu et all 
(2010)). These priors are wide enough to ensure that CMB con- 
straints are not double counted when our result s are combined with 
CMB data dChuangl Wan g. & Hemanthdl2012l) . We also marginal- 
ize over the amplitude of the galaxy correlation function, effec- 
tively marginalizing over a linear galaxy bias. 



4 RESULTS 

4.1 Measurement of multipoles 

FigUHUH] and [4] show the effective monopole (£o), quadrupole 
(£2), hexadecapole (£4), and hexacontatetrapole (£e) of measured 
from SDSS LRGs compared with the averaged effetive multipoles 
measure d from the mock catalo gs. We are using the same scale 
range as lChuang & WangH201 lh (s = 40-120 ft _1 Mpc) for com- 
parison and the bin size used is 5 /i _1 Mpc. The data points from 
the multipoles in the scale range considered are combined to form 
a vector, X (see equationl[T9ll). 

We find that £t(s) deviates on scales of s < 60Mpc//i from 
the measurement from mock data (in contrast to £o(s), &(s), and 
leading to a significant difference in the measured mean 
values of H(z), Da(z), and f3(z) from £0 + £2 ar, d £0 + £2 + £4, 
thus it should not be used in deriving parameter constraints. 

A frequently used combination of the monopole and the 
quadupole is the normalized quadrupole, defined by 



Us) 



£o(s)-(3/ S 3 ) j;Us')s' 2 ds'' 



(23) 



For comparison with previous work, we measure the effective 
normalized quadrupole defined by 



Q(s) 



Us) 



(24) 



|o( S )-(3/ S 3) J2 Us')s' 2 As 

0<s'^s 

from SDSS LRGs and show it in Fig(5] It is in good agreement with 
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Figure 2. Measurement of quadrupole of the correlation function of 
the SDSS DR7 LRG (diamond data points) comparing to the averaged 
quadrupole of the correlation functions of the mock catalogs (solid line). 
The error bars are taken as the square roots of the diagnal elements of the 
covariance matrix. 



Figure 4. Measurement of hexacontatetrapole of the correlation function 
of the SDSS DR7 LRG (diamond data points) comparing to the averaged 
hexacontatetrapole of the correlation functions of the mock catalogs (solid 
line). The error bars are taken as the square roots of the diagnal elements of 
the covariance matrix. 




so 

s (Mpc/h) 



Figure 3. Measurement of hexadecapole of the correlation function of the 
SDSS DR7 LRG (diamond data points) comparing to the averaged hexade- 
capole of the correlation functions of the mock catalogs (solid line). The 
error bars are taken as the square roots of the diagnal elements of the co- 
variance matrix. 



the expe ctation from the LasDam as mocks, as well as with previous 
work by Samush ia et al. I d201ll) . 




s (Mpc/h) 



Figure 5. Measurement of the normalized quadrupole from the SDSS DR7 
LRG (diamond data points) comparing to the mean measurement from the 
mock catalogs (solid line). The error bars are taken as the square roots of 
the diagnal elements of the covariance matrix. Dashed line is the theoretical 
prediction for beta = 0.3 assuming linear power spectrum and small-angle 
approximation. 



4.2 Measurement of {#(0.35), D A (0.35), /3(0.35)} 

We now present the model independent measurements of 
the parameters {#(0.35), Da(0.35), Q m h 2 ,f3}, obtained by 
using the method described in previous sections. We also 
present the derived parameters including #(0.35) r s {zd)/c and 
Da (0.35) /r s (z,i) which have tigher constraints. 

First, we apply our method of using multipoles of the corre- 
lation function to the individual LasDamas mock catalogs and find 
that our measurements are consistent with the input parameters of 



the simulations and also the results from using full 2D 2PCF. We 
also find that the constraints from using monopole + quadrupole + 
hexadecapole is almost as good as using full 2D correlation func- 
tion, but using extra hexacontatetrapole does not provide significant 
improvement to the constraints. Then, we apply our method on the 
observational data. Finally, we compare our results with other pre- 
vious work by restricting our method to the same conditions ac- 
cordingly. 
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mean 


(7 


input value 


H(0 35) 


81.1 


5.6 


81.79 


D4(0 35) 


1017 


63 


1032.8 




0.119 


0.014 


0.1225 





0.325 


0.076 




H(0.35)r s (z d )/c 


0.0436 


0.0030 


0.0434 


D A (0.35)/r a (z d ) 


6.29 


0.36 


6.48 



Table 1. The mean and standard deviation of the distribution of the mea- 
sured values of {if(0.35), D A (0.35), n m h 2 , 0, H(0.35)r s (z d )/c, 
Da{0.35) /r s (z d )} from each monopole + quadrupole (£o + £2) of 40 
LasDamas mock catalogs (which are indexed with 01a-40a). Our measure- 
ments are consistent with the input values within lcr, where each a is com- 
puted from the 40 means measured from the 40 mock catalogs. The unit of 
# is km s _1 Mpc -1 . The unit of Da and rs(Zd) is Mpc. 



4.2.1 Validation Using Mock Catalogs 

In order to validate our method, we have applied it to the first 
40 LasDamas mock catalogs (which are indexed with 01a-40ajfl 
Again, we apply the flat and wide priors (±7awMAP7) on Qbh 2 
and n s , centered on the input values of the simulation (Qth 2 = 
0.0196 and n 3 = 1). 

Table Q] shows the means and standard deviations of the dis- 
tributions of our measurements of {#(0.35), Z?a(0.35), Q m h 2 , 
P, H (0.35) r s {zd)/c, DA(0.35)/r a (zd)} from each monopole 
+ quadrupole (fn + £2) of the LasDamas mock catalogs of the 
SDSS LRG sample. Table [2] and [3] show the measurements us- 
ing monopole + quadrupole + hexadecapole (£0 + £2 + £4) and 
monopole + quadrupole + hexadecapole + hexacontatetrapole (£0 + 
£2 + £4 + £e) respectively. For c omparison, Table [4] show s the re- 
sults using full 2D method from IChuang & Wand ( 1201 lh . but we 
are using only 40 for this study instead of 80 mock catalogs. These 
are consistent with the input parameters, establishing the validity 
of our method. 

While the constraints from using £0 +£2 +£4 are significantly 
tighter than using £0 + £2, the constraints from using £0 + £2 + 
£4 + £e are nearly the same as that from using £0 + £2 + £4- This 
indicates that £0 + £2 +£4 captures nearly all of the information that 
can be extracted from the data given the noise level. In principle, 
one could obtain better constraints by including more multipoles. 
However, the tradeoff is introduing noise to the covariance matrix 
which could be a problem, since the number of the mock catalogs 
used to construct the covariance matrix is not big enough. We also 
show the measurements of #(0.35) r s (z d )/c, Da (0.35) /r 3 (z d ), 
Q m h 2 , and /? of each mock catalog in Fig. [6] Fig. [7] Fig. [8] and 
Fig.|9]to show the scattering among different mock catalogs and the 
deviations among different methods. One can see the measurements 
among different methods are consistent for most mock catalogs, but 
there are still some obvious deviations (> la) for a few cases. 

An important point to note is that since the mock data do not 
include unknown systematic effects, the mean values of estimated 
parameters remain nearly unchanged as more multipoles measured 
from data are added to the analysis and the parameter constraints 
are tightened with the addition of information. 



3 We only use 40 instead of 160 mock catalogs because the MCMC is com- 
putationlly expensive. However, the covariance matrix is contracted with 
160 mock catalogs. 





mean 


(7 


input value 


# (0 35) 


80.4 


4.9 


81.79 


Da(0 35) 


1027 


56 


1032.8 




0.116 


0.013 


0.1225 


P 


0.327 


0.066 




H(0.3S)r 3 (z d )/c 


0.0435 


0.0025 


0.0434 


D A (0.35)/r a (z d ) 


6.31 


0.31 


6.48 



Table 2. The mean and standard deviation of the distribution of the mea- 
sured values of {#(0.35), D A (0.35), Q m h 2 , 0, #(0.35) r a (z d )/c, 
D A (0.35) /r a (z d )} from each monopole + quadrupole + hexadecapole 
(£0 + £2 + £4) of 40 LasDamas mock catalogs (which are indexed with 
01a-40a). Our measurements are consistent with the input values within lcr, 
where each a is computed from the 40 means measured from the 40 mock 
catalogs. The unit of # is km s — 1 Mpc . The unit of Da and r a (z d ) is 
Mpc. 





mean 


(7 


input value 


#(0.35) 


80.3 


5.1 


81.79 


D A (0.35) 


1021 


48 


1032.8 




0.116 


0.013 


0.1225 





0.324 


0.075 




H(0.35)r a (z d )/c 


0.0434 


0.0026 


0.0434 


D A (0.35)/r a (z d ) 


6.28 


0.26 


6.48 



Table 3. The mean and standard deviation of the distribution of the mea- 
sured values of {#(0.35), £> A (0.35), Q m h 2 , 0, #(0.35) r a (z d )/c, 
Da(0. 35) /r a (z d )} from each monopole + quadrupole + hexadecapole + 
hexacontatetrapole (go + £2 + £4 + of 40 LasDamas mock catalogs 
(which are indexed with 01a-40a). Our measurements are consistent with 
the input values within lcr, where each cr is computed from the 40 means 
measured from the 40 mock catalogs. The unit of # is kms" 1 Mpc . 
The unit of Da and r a (z d ) is Mpc. 

4.2.2 Measurements from SDSS DR7 LRG 

Table [5] lists the mean, rms variance, and 68% confidence 
level limits of the parameters, {#"(0.35), Da(0.35), Q m h 2 , ft, 
#"(0.35) r 3 {z d )/c, Da(0.35) /r s (z d )}, derived in an MCMC like- 
lihood analysis from the measured £0 + 62 of the correlation func- 
tion of the SDSS LRG sample. Table[6]lists the mean, rms variance, 





mean 


a 


input value 


#(0.35) 


80.5 


4.0 


81.79 


D A (0.35) 


1026 


54 


1032.8 




0.117 


0.013 


0.1225 





0.353 


0.060 




#(0.35) r B (z d )/c 


0.0435 


0.0021 


0.0434 


D A (0.35)/r s (z d ) 


6.31 


0.26 


6.48 



Table 4. The mean and standard deviation of the distribution of the mea- 
sured values of {#(0.35), D A (0.35), Q m h 2 , 0, #(0.35) r s (z d )/c, 
Da(0.35) /r s (z d )} from each 2D 2PCF of 40 LasDamas mock catalogs 
(which are indexed with 01a-40a). Our measurements are consistent with 
the input values within lcr, where each a is computed from the 40 means 
measured from the 40 mock catalogs. The unit of # is kms -1 Mpc -1 . 
The unit of Da and r s (z d ) is Mpc. 
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Figure 6. Measurements of the means and standard deviation of 
H(0.35) r s (z ( i)/c from 40 individual mock catalogs (indexed as Ola to 
40a). The blue circles show the measurements using §o + ?2- The red din- 
monds show the measurements using £o ± £2 + £4 ■ The black crosses show 
the measurements using full 2D 2PCF method from our previous work. The 
black line shows the theoretical value computed with the input parameters 
of the simulations. 



Figure 8. Measurements of the means of f2 m b? from 40 individual mock 
catalogs (indexed as Ola to 40a). The blue circles show the measurements 
using £o ±?2 ■ The red dinmonds show the measurements using £o +£2 +£4 ■ 
The black crosses show the measurements using full 2D 2PCF method from 
our previous work. The black line shows the theoretical value computed 
with the input parameters of the simulations. 
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Figure 7. Measurements of the means of DA(0-35)r s (zd) from 40 indi- 
vidual mock catalogs (indexed as Ola to 40a). The blue circles show the 
measurements using £o + £2 ■ The red dinmonds show the measurements 
using £o + £2 + ?4- The black crosses show the measurements using full 
2D 2PCF method from our previous work. The black line shows the theo- 
retical value computed with the input parameters of the simulations. 

and 68% confidence level limits of the same parameter set from the 
measured £0 + £2 + £4 of the correlation function of the SDSS 
LRG sample for this parameter set. The \ P er degree of freedom 
(X 2 / d -°-f-) is 1.23 for | + £2 and is 1.06 for £ + £ 2 + £4. These 
are independent of a dark energy model, and obtained without as- 
suming a flat Universe. The measuements from these two methods 
have deviations which are close to lcr which occur for some mock 
catalogs as well (i.e. the result from the 12th mock catalog in Fig. 
[6}. The low measurement of Q. m h 2 and high measurement of /3 are 
found from the mock catalogs, too (i.e. result from the 16th mock 
catalog in Fig.[8]and Fig- [9j - While these cases are rare among the 



Figure 9. Measurements of the means of ft from 40 individual mock cata- 
logs (indexed as Ola to 40a). The blue circles show the measurements using 
f + £2 ■ The red dinmonds show the measurements using £0 + £2 + £4 ■ The 
black crosses show the measurements using full 2D 2PCF method from our 
previous work. 

mock catalogs, it is possible that the measurements from the ob- 
servational data are distorted by some unknown systematic error 
and not just sample variance. These deviations might provide not 
only the hints of the appearance of the systematic errors but also 
the clues for studying them. Here, to be conservative, we use the 
measurements using £0 + £2 as our fiducial results, since £4 has 
obvious deviation from the measurement from the mock data (in 
contrast to £0, £2, and see Figs. 1-4). 

Table UJ gives the normalized covariance matrix for this pa- 
rameter set measured using £0 + £2- While the measurement of j3, 
0.44 ±0.15, is higher than what we expect (i.e/3 = 0.325 ±0.076 
from the mock catalogs using £0 + £2), note that there is a negative 
correlation between /3 and Q m h 2 and the correlation coefficient 



© 0000 RAS, MNRAS 000, 000-000 



8 Chuang & Wang 





mean 


<7 


lower 


upper 


HYO 3^1 


79.6 


8.8 


70.9 


87.8 


D A (0.35) 


1060 


92 


970 


1150 




0.103 


0.015 


0.088 


0.118 





0.44 


0.15 


0.29 


0.59 


#(0.35) r s (z d )/c 


0.0435 


0.0045 


0.0391 


0.0477 


D A (0.35)/r s (z d ) 


6.44 


0.51 


5.99 


6.90 



Table 5. The mean, standard deviation, and the 68% C.L. 
bounds of {#(0.35), D A (0.35), n m h 2 , 0, #(0.35) r s (z d )/c, 
D A (0.35)/r 3 (z d )} from SDSS DR7 LRGs using £ + 6- The unit of # 
is kms~ 1 Mpc - . The unit of D A and r a (z d ) is Mpc. 





mean 


<7 


lower 


upper 


#(0.35) 


87.3 


6.7 


80.8 


93.7 


D A (0.35) 


1095 


59 


1037 


1153 


^ m h 


0.107 


0.015 


0.093 


0.122 





0.54 


0.11 


0.44 


0.65 


H (0.35) r a (z d )/c 


0.0472 


0.0033 


0.0441 


0.0503 


D A (0.35)/r s (z d ) 


6.75 


0.25 


6.52 


6.98 



? 1 




D A /D A( . d 

Figure 11. 2D marginalized contours (95% C.L. ) for fl^tz)/!)':" and 
H(z)/H true for the comparison with fig. 6 in IKazin et ail 120121) . The 
black solid contour is measured using £o + £2 and the blue dotted contour 
is measured using gp + £2 + £4 ■ Our constraints are similar with the results 
in lKazin et all l2012h . 



Table 6. The mean, standard deviation, and the 68% C.L. 
bounds of {#(0.35), D A (0.35), Q m h 2 , 0, H (0. 35) r s (z d )/c, 
D A (0.35)/r s (z d )} from SDSS DR7 LRGs using £0 +6 +I4. The unit 
of # is km s — 1 Mpc -1 . The unit of D A and r s (z d ) is Mpc. 

is -0.2549. Thus the somewhat high j3 value is mildly correlated 
with the somewhat low Q, m h 2 value. The most robust measure- 
ment s are that of \ H (0.35) r s (z d )/c, D A (0.35)/r s (z d )}, same 
as in IChuang&Wan3 l l201ll) . These can be used to combine with 
other data se ts and constraining dark energy and cosmological pa- 
rameters, see Wang, Chuang & Mukherjee (2012). 

Fig. [TO] shows one and two-dimensional marginalized con- 
tours of the parameters, {#(0.35), Da (0.35), Q. m h 2 , (3, 
H(0.35)r 3 (z d )/c, D A (0.35)/r 3 (z d )}, derived in an MCMC like- 
lihood analysis from the measured £0 + £2 of the SDSS LRG sam- 
ple. 

4.3 Comparison with Previous Work 

While we have developed a general method to measure the dark 
energy and cosmological parameters which could be extracted from 
the galaxy clustering data alone, we restrict our method now by 
fixing some parameters to obtain the results for comparison with 
the pr evious work by o thers. 

IKazin et alj (20121) present some tests for measuring # ( z ) and 
D A (z) using averaged multipoles of the correlation function from 
the LasDamas mock catalogs. They assume a larger survey volume 
(~ 12 times) by dividing the covariance matrix by \/160. They 
use the averaged mulitpoles of the correlation function from the 
mock catalogs as the theoretical model, which is equivalent to fix- 
ing Q, m h 2 , flbh 2 , and n s to the input parameters of the simulations. 
We fix the damping factor, k+ — 0.13 /iMpc -1 , and the pairwise 
peculiar velocity dispersion, a v = 300 s -1 km, which give a good 
fit to the averaged correlation function of t he mock catalogs. Cor- 
responding to the bottom panel of fig. 6 in.Kaz in et alj d2012h . we 
measure the hubble parameter and angular diameter distance by 
marginalizing over the amplitude of the correlation function and 



the linear redshift distortion parameter and using the scale range, 
s = 40 — 150/i -1 Mpc. Our results are shown in Fig. Illlwhich is 
similar to theirs. The 1-D marginalized uncertainties of {#, D A , 
P} we measure using £ +6 +£4 is {1.1 7%, 0.81%, 4.45%} which 
are similar to their resu lts, {1.42% 0.76% , 4.95%} (the numbers 
are taken from Fig. 7 in lKazinetal.U2012h l They derive the theo- 
retical multipoles analytically, instead of using the same definition 
applied to the observational data. In principle, it could introduce bi- 
ases to the measurments. However, the effect might be minimized 
since they contract the theoretical model based on the measured 
multipoles from the mock catalogs that is equivalent to computing 
the theoretical multipoles with the same definition applied to the 
obse rvational data. 

ICabre & Gaztanagal d2009l) measure /3 from SDSS DR6 LRG 
using the normalized quadrupole defined by Eq. d23) , To compare 
with their results, we make similar assumptions, and use monopole- 
quadrupole method with fixing Q m = 0.25, Qb = 0.045, h = 
0.72, n s = 0.98, fc* = 0.11, and a v = 300km/s in the ACDM 
model (#(0.35) and Da(0.35) would also be fixed accordingly). 
Considering the scale range, s = 40 - lOO/i^Mpc, we obtain 
ft — 0.333 ± 0.055, in excellent agreement with their meausure- 
ment of /3 = 0.34 ± 0.05. Since the definition of the normal- 
ized quadrupole includes a integral of monopole with the minimum 
boundary from s = 0, the advantage of using our effective multi- 
pole method instead of normalized quadrupole method is to avoid 
the distortion from the small scales where the scale dependent un- 
certainties are not well known. However, the distortion might be 
negligible compared to the statistical uncertainty of current mea- 
surements^ 

ISong et al. I l l201lh split the same galaxy sample (SDSS DR7 
LRG) to two redshift slides and obtained f3(z — 0.25) = 
0.30t^ols and P( z = °- 38 ) = °- 39 ± °- 056 without consider- 
ing the geometric distortion. Since we use the whole galaxy sample 
as one redshift slice, our results are not directly comparable with 
theirs. 
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Figure 10. 2D marginalized contours (68% and 95% C.L.) for {#(0.35), D A (0.35), tt m h 2 , ff(0.35) r s (z d )/c, D A (0.35)/r a (z d )}. The diagonal 
panels represent the marginalized probabilities. The unit of H is kms -1 Mpc -1 . The unit of Da and r s (zd) is Mpc. 



5 SYSTEMATIC TESTS 

Table. [8]shows the systematic tests that we have carried out varying 
key assumptions made in our analysis. These include the multipoles 
used, the range of scales used, and the bin size used to calculate the 
correlation function. 

We use the results using £n + £2 as our fiducial results. We 
find that the constraints are stronger while using £0 + £2 + £4, but 
using £0 + £2 + does not improve the constraints siginicantly. 
Therefore, it seems that £n + £2 + £4 contains most of the informa- 



tion from the 2D 2PCF. However, since we suspect there is some 
systematic error for the measurement of £4 from observational data 
(see Fig. 3), we do not include it for our fiducial results. 

We vary the scale chosen and the bin sized used and find the 
results are basically consistent. However, we find the measurement 
of DA{z)/r s (zd) is more stable than that of H(0.35) r 3 (zd)/c. 
It might indicate the appearance of the systematic error from the 
measurement of the correlation function in the direction along the 
line of sight. 



© 0000 RAS, MNRAS 000, 000-000 



10 Chuang & Wang 





ff(0.35) 


D4(0.35) 







H(0.35)r B (z d )/c 


D A (0.35)/r s (z d ) 


ff(0.35) 


1 


0.2669 


0.3529 


0.5802 


0.9259 


0.4832 


-0.4(0.35) 


0.2669 


1 


-0.3835 


0.6307 


0.4586 


0.8814 




0.3529 


-0.3835 


1 


-0.2549 


-0.0057 


0.07 


/8 


0.5802 


0.6307 


-0.2549 


1 


0.7176 


0.575 


ff(0.35)r s ( Z<i )/c 


0.9259 


0.4586 


-0.0057 


0.7176 


1 


0.4981 


D A (0.35)/r s (z d ) 


0.4832 


0.8814 


0.07 


0.575 


0.4981 


1 



Table 7. Normalized covariance matrix of the measured and derived parameters, {.ff (0.35), 1)^(0.35), Cl m h 2 , /3, H(0.35) r s (z l i)/c, Z}a(0.35)/V s (z d )} 
from SDSS DR7 LRGs using £0 + 6. 



fl"(0.35) D A (0.35) 



f) H(0.35)r s (z d )/c D A (0.35)/r s (z d ) 



£o + €2 (fiducial) 
f + €2 + £4 

fo + €2 + €e 
30 < s < 120 
50 < s < 120 
40 < s < 110 
40 < s < 130 
bin size = 4 h ~ 1 Mpc 
bin size = 8 h ~ 1 Mpc 



n+8.2 
-8.6 
87.3 ± 6.4 
r + 8.7 
J -8.9 
85. 1 4 



78., 

-8.2 
77 5+ 8 - 2 

73 5 + 6 7 

83.5+|-| 
+ 7.'l 

-7.2 
+ 8.1 
8.6 



77.3 
79.1 



1063+|J 
1095 ± 58 

1025112 
1072+^ 

1034+ 103 



1064 4 



-109 
^ + 77 
-76 

1082+1/? 
— 75 

1026+1® 

1 ncr+84 

1065 -78 



0.103 ± 0.015 

107+ ' 015 
-0.014 

0.107 ± 0.016 

0.115 ± 0.014 

101+ ' 019 
" " -0 .018 

0.107 ± 0.015 
0.105 ± 0.015 
0.113 ± 0.017 
0.112 ± 0.017 



45+ ' 15 
° -0.14 

0.54 ± 0.11 

0.41 ± 0.14 

0.38 ± 0.10 

5+ 019 
u - 5 -0.18 

0.35 ± 0.11 
0.48 ± 0.15 
0.39 ± 0.12 
0.42 ± 0.13 



0437+ 0041 
" -0.OO43 

0.0472 ± 0.0031 

0424+ 0044 
-0.0043 

0453+ 0037 
U.U4&3_ 0039 

0.0425±°;°°» 

0.0398+.^ 

0454+ 0041 
UUJ -0.0043 

0.0413 ± 0.0034 

0423+ - 0039 
u - u4 ^-0.0041 



6.48 1 



a+0.44 

-0.43 
6 75+ ' 24 

6.31 ± 0.49 

6 71+ 031 
-0.30 

6 27+ 55 
"'^'-0.61 

6 54+ ' 41 
u - J *-0.40 

6 63+ ' 37 
6.39 ± 0.45 



Table 8. This table shows the systematic tests with using diffirent combination of multipoles, the scale range, and the bin size used. The fiducial results are 
obtained by considering £o + £2 + £4, the scale range (40 < s < 120 h~ 1 Mpc), and the bin size (5h~ 1 Mpc). The other results are calculated with only 
specified quantities different from the fiducial one. The unit of H is kms -1 Mpc -1 . The unit of Da and r a (z l i) is Mpc. 



6 CONCLUSION AND DISCUSSION 

We have demonstrated the feasibility of using multipoles of the cor- 
relation function to measure H(z), Da(z), Q, m h 2 , and /3 by apply- 
ing the method to the individual mock catalogs from LasDamas in 
an MCMC likelihood analysis. 

The method we developed is modified from Chuang & Wang 

1) which is the first application of including the geo- 
matric distortion (also kn own as Alcock-Paczynski test, see 
lAlcock & Paczvnskl d 19791) ) on the galaxy clustering data at large 
scales. We compute the multipoles from the theoretical and ob- 
served 2D 2PCF in the same way, thus no more approximation is 
made besides the distance of any pair of galaxies can be converted 
with two stretch factors between different models in the redshift 
range considered. 

We have also obtained the constraints for the mea- 
sured and derived parameters, {ff(0.35), Da(0.35), Q. m h 2 , ft, 
.ff (0.35) r a (zd)/c, Z)A(0.35)/r s (zd)}, from the multipoles of the 
correlation function from the sample of SDSS DR7 LRGs which 
are summarized by Tables |5]and|6] 

We find that while the mean values of estimated parame- 
ters remain stable (with rare deviations) for the mock data when 
higher multipoles are used, this is not true for the SDSS DR7 LRG 
data. We find #(0.35) r s (z d )/c = 0.0437tg ^ using monopole 
+ quadrupole, and #(0.35) r a {z d )/c = 0.0472 ± 0.0031 using 
monopole + quadrupole + hexadecapole. This deviation is caused 
by the hexadecapole measured from SDSS LRGs deviating from 
the value expected from the mocks on scales less than about 
60Mpc//i. This suggests that the multipole method can be used 
to isolate sytematic uncertainties in the data, and provide a useful 
cross-check of parameter measurements from the full correlation 
function. 
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